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The properties of two forms of the gradient expanded Kadanoff-Baym equations, 
i.e. the Kadanoff-Baym and Botermans-Malfliet forms, suitable to describe the 
transport dynamics of particles and resonances with broad spectral widths, are dis- 
. cussed in context of conservation laws, the definition of a kinetic entropy and the 

o , 

, possibility of numerical realization. Recent results on exact conservations of charge 

' and energy-momentum within Kadanoff-Baym form of quantum kinetics based on 

local coupling schemes are extended to two cases relevant in many applications. 
These concern the interaction via a finite range potential, and, relevant in nuclear 
I and hadron physics, e.g. for the pion-nucleon interaction, the case of derivative 
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coupling. 



INTRODUCTION 



CJ ' Ever since L. Boltzmann had suggested his famous kinetic equation, the field of non- 

^ ■ equihbrium physics and stochastic processes has grown tremendously, expanding into various 

. ^ , directions. The interactions among particles driven by mean fields were included, quasipar- 
' tides were introduced in order to include much of the medium effects, the kinematics was 
extended to the relativistic case, ultimately theoretical foundations of transport equation 
were given from an underlying quantum many-body or field theory. In this line of achieve- 
ments also stands the work of G.I. Budker and S.T. Belyaev who demonstrated the Lorentz 
invariance of the relativistic distribution function and derived relativistic Fokker-Planck ki- 
netic equation The work entered into many textbooks and found numerous applications 
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in atomic physics and electron-positron plasma. Presently the relativistic transport concepts 
are a conventional tool to analyze the dynamics of dense and highly excited matter produced 
in relativistic heavy-ion collisions. 

Along with semiphenomenological extensions a great progress was achieved in microscopic 
foundation of the kinetic theory, which is mainly associated with the names of Bogolyubov, 
Born, Green, Kirkwood, Yvon, and Zubarev. The appropriate frame for description of non- 
equilibrium processes within the real-time formalism of quantum-field theory was developed 
by Schwinger, Kadanoff, Baym and Keldysh 0,0,0]. This formalism allows extensions of 
the kinetic picture beyond conventional approximations (like the quasiparticle one) and has 
found now numerous applications in many domains of physics. 

The interest in transport descriptions of heavy-ion collisions beyond the quasiparticle 
approximation was initiated by Danielewicz using the gradient expanded Kadanoff- 



it ejcpanded Kadanoii- 



Baym equations. These attempts have recently has been revived 

Q] in order to properly describe the transport properties of broad resonances (like the p 
meson and A isobar). In the dense environment also stable particle acquire a considerable 
width because of collisional broadening. A proper dynamical treatment of their widths 
in the dense nuclear medium within a transport theoretical concepts is still a challenging 
problem. Transport approaches for treating such off-shell dynamics were proposed in refs. 

. They all were based on the Kadanoff-Baym equations which 
describe the non-equilibrium quantum evolution at the truncation level of the Schwinger- 
Dyson equation. Expanded up to the first space-time gradients these provide equations 
for the one-body phase-space distribution functions with collision term and Poisson-bracket 
terms arising from the first order gradient terms. Presently two slightly different forms of the 
gradient-expanded Kadanoff-Baym equations are used: the original Kadanoff-Baym (KB) 
form |3;] , as it follows right after the gradient expansion without any further approximations, 
and the Botermans-Malfliet (BM) one , which is derived from the KB form by omitting 
certain second-order space-time gradient corrections. 

In this paper we would like to compare these two forms of "quantum" kinetic equations 
and discuss their advantages and disadvantages from the point of view of their conserving 
properties, the possibility of numerical realization, etc., sect. |n] Technical details on the 
conserving-properties are deferred to the appendices (jBHD|) . since they illustrate some of the 
general consideration of ref. together with some extensions to cases particularly relevant 
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in nuclear physics. Append. O treats nonrelativistic nucleon-nucleon interactions via a 
potential of finite range. The derivative coupling is considered in Append. |D]at the example 
of the P-wave pion-nucleon interaction. In sect. IIIH we supplement some considerations 
about the construction of a kinetic entropy within these two transport schemes. To make 
the paper self-contained, we summarize the time-contour matrix notation in Append. 
and introduce the ^-functional formalism for derivative coupling in Append. El A summary 
is given in sect. IIVI 



II. OFF-SHELL KB AND BM KINETICS 



In this section we summarize the formulation of the off-shell kinetic equations in the two 
different forms: in the KB form, i.e. as it follows right after the gradient expansion of the 
exact Kadanoff-Baym equations, and in the BM form which differs from the KB form 
only in the second order of the gradient expansion. We assume the reader to be familiar 
with the real-time formulation of non-equilibrium many-body theory and use the contour 
matrix notation, detailed in Append. El 

Starting point of all considerations is the set of Kadanoff-Baym equations which express 
the space-time changes of the Wigner transformed correlation function iG^^{X,p) in terms 
of the real-time contour convolution of the self-energy S with the Green function G. We 
give the kinetic equation in compact notation (cf. Eq. ()A4j) ) 

v,d'^iG-+{X,p) = [J:^G-G®^-^l with v^ = -^Go\p), (1) 

where Gq^{p) is the Fourier transform of the inverse free Green function 

, — fTi^ for relativistic bosons 

Go\p) = { ^ (2) 

I Po ~ p i\2m) for non-rel. fermions or bosons. 

For a complete definition, Eq. has to be supplemented with further equations, e.g., for 
the retarded Green function together with the retarded relations ()A6|) . In addition to these 
equations the exact set of Kadanoff-Baym equations also includes the prototype of the mass- 
shell equation, which we also discuss below. If a system under consideration is only slightly 
spatially inhomogeneous and slowly evolving in time, a good approximation is provided by 
an expansion up to first order in space-time gradients. Then the main problem to arrive at 
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a proper kinetic equation consists in accurately disentangling the rather complicated r.h.s. 
of Eq. (0). This problem in the context of conserving approximations will be addressed 
here. 



A. <l>-derivable approximations 

In actual calculations one has to use certain approximations or truncation schemes to 
the exact non-equilibrium theory, which make conserving properties (such as charge and 
energy-momentum conservations) and thermodynamic consistency of the transport theory 



not evident. It was shown 
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17| that there exist a class of self-consistent approximations, 
called $ derivable approximations, which are conserving at the expectation value level and 
at the same time thermodynamically consistent, i.e. they provide true Noether currents and 
a conserved energy-momentum tensor. In these schemes the self-consistent self-energies are 
generated from a functional ^[G] through the following variational procedure 

'E (X ) T ^^'^^^^ I ^ ^^^^ ^^^^^ ■ ke{ +} (3) 

' 5iG"=^(X,p) |i for complex fields ' 

The functional specifies the truncation scheme. It consists of a set of properly chosen 
closed two-particle irreducible diagrams, where lines denote the self-consistent propagator 
G, while vertices are bare. The functional variation with respect to G diagrammatically 
implies an opening of a propagator line of 

Particular examples of $ derivable approximations can be found in Appendices IHl and \D\ 
which consider applications of the general formalism to cases important in nuclear physics. 
The treatment of the pion-nucleon derivative coupling in Append. O requires the corre- 
sponding extension of the $ derivable formalism, which has not been done up to now. There- 
fore, in Append. El we perform such an extension and derive the relevant modifications of 
the variational rules and the ensuing additional terms in the current and energy-momentum 
tensor expressions. 

The conserving properties of these approximations are exact at the level of Kadanoff- 
Baym equations (^, while after the expansion up to the first space-time gradients they are 
generally expected to be only approximately fulfilled. 
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B. Physical notation 

It is helpful to eliminate the imaginary factors inherent in the standard Green function 
formulation and introduce quantities which are real and, in the quasi-homogeneous limit, 
positive, with clear physical meaning, thereby. Thus instead of Green functions G''^{X,p) 
and self-energies p) with i,j G { — h} (see Append. El in the Wigner representation 

we use the kinetic notation of ref. We define^ 

F{X,p) = A{X,p)f{X,p) = {T)iG-^{X,p), 

F{X,p) = A{X,p)[l T fiX,p)] = iG+-(X,p), (4) 

for the generalized Wigner functions F and F and the corresponding /our-phase-space dis- 
tribution functions f{X,p) and Fermi/Bose factors [1 =F f{X,p)]. Here 

A{X, p) = -2Im G^{X, p) = F±F (5) 

is the spectral function, and is the retarded propagator. The spectral function satisfies 
the sum rule 

[ — ^A(X, = 1 for nonrelativ. particles, (6) 
Jo 27r 

/--^PqA{X,p) = 1 for relativ. bosons, (7) 

which follows from the canonical equal-time (anti) commutation relations for (fermionic) 
bosonic field operators. Likewise the gain and loss rates of the collision integral are defined 

as 

r-{x,p) = Ti^-+{x,p), r-^(x,p) = is+-(x,p) (8) 

with the damping width 

T{X,p) = -2ImS^(X,p) = T°''\X,p) ± r"^(X,p), (9) 

where is the retarded self-energy. 

In terms of above kinetic notation, the gradient-expanded Kadanoff-Baym equations are 
reduced to equations for four real quantities: two equations for the real and imaginary parts 



^ Here and below the upper sign corresponds to fermions, while the lower sign, to bosons. 
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of the retarded Green function, while there are two equations for the phase-space occupation 
F: the KB kinetic equation and the prototype "mass-shell equation" . The latter doubling 
of equations reflects the well known redundancy of the Kadanoff-Baym equations. Before 
the gradient expansion both equations are completely identical. However, after the gradient 
expansion their interrelation is no longer obvious and deserves special care (see below). 

The equations for the retarded pro pag ator in first-order gradient approximation can be 
immediately solved with the result 0, llS^ 

qR 

Mix r>\ ^^\'^ \ n W V 

ReG'^(X,p) 



M2(X,p) + r2(X,p)/4 
M(X,p) 



M2(X,p) + r2(X,p)/4 



(10) 



with the "mass" function 



M(X,p) = Go^(p)-ReE«(X,p). 



(11) 



Although the solution (fTU|) is simply algebraic, it is vahd up to first-order gradients. 



C. Kadanoff Baym form 

In terms of above notation, the KB kinetic equation for F in the first-order gradient 
approximation takes the form 

VF{X,p) - {r-,ReG^} = C{X,p). (12) 

We denote this as the quantum transport equation in the KB-choice^ . Here the differential 
drift operator is defined as 

and {...,...} denotes the four-dimensional Poisson bracket 



^ If the system consists of several different particle species, there is a set of coupled kinetic equations 
corresponding to each species (e.g., see Append. ID|| . 
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in covariant notation. Please note that now after gradient approximation all quantities on 
the l.h.s. are to be taken in local approximation, i.e. void of any further gradient terms. 
Thus, the occurring self-energies are obtained evaluating the diagrams as in momentum 
representation with the coordinates X of all Greens functions kept identical. The r.h.s. 
specifies the collision term^ 

c{x,p) = r-{x,p)F{x,p) - r^-\x,p)F{x,p). (15) 

If the diagrams for the self-energy contain internal vertices, which give rise to memory or 
non-local effects, the gain and loss rates contain additional gradient terms which have to be 
constructed, e.g., according to the rules given in [l^. The resulting local part of the collision 
term is charge (e.g., the baryonic number) and energy-momentum conserving by itself 

= 0. (16) 

Here and below, e denotes the elementary charge, while Tr implies the sum over all possible 
internal degrees of freedom, like spin, and over possible particle species. We do not explicitly 
introduce the particle-species label to avoid overcomplication of equations. In terms of a 
local functional the explicit form of the local collision term is 

cf. Eq. (jHJ. In this paper we limit the considerations to cases void of memory effects in this 
collision term. The latter effects were studied in 
Relation (fTB|) permits us to derive the current 

d^p 




j^B-esiX) = eTr J 



a TP r)F" 
v^'FiX, p) + Re S^-^ - Re 

OPu OPn 



(18) 



of a charge e (e.g., the baryonic one) from the KB kinetic equation (fT^ which is conserved 

«B-efr(^) = 0. (19) 
Note that this current formally differs from the true Noether current 

= eTr I ^,v^F{X,p)+j^^^^^{X), (20) 



See an example in Eq. (|C3I 
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which follows right from the operator expression for this quantity, cf. ref. [8| and Appendix 
IbI The additional term j'^^^^^ appears only in the case of derivative couphng, see Eq. ()B19|) . 
In view of the gradient approximation employed, one could generally expect both currents to 
differ beyond the validity range of the gradient approximation. However, as demonstrated 
in detail in ref. these two currents are exactly equal for ^-derivable approximations, 
if a consistent gradient expansion is performed also in the gain and loss rates (jH)) of the 
collision term (1151) . In this case the exact conservation of the Noether current results from 

n 

the corresponding invariance of the $ functional - (Eq. (6.9) in JJ}) ^ which survives the 
gradient expansion 

J (^ KReS^,F} - {ReG^,r-} +C] = d,ji^,^^y (21) 

The latter relation is written down for the general case of memory or non-local effects 
included in C. If such effects in C are absent, the collision term drops out of Eq. (pT|l 
according to Eq. (fT^ . 

Within ^-derivable approximations also the conservation of energy-momentum can be 
established, for local (point-like) couplings providing a local energy-momentum tensor. The 

weighted four-momentum integral of the KB kinetic equation leads to the following con- 
sistency relation j^"^ 

= (^).., / 1^ [{ReS«,F} - {ReG-,T-} + C] , (22) 

which is again exact after the gradient expansion as shown in ref. Q| (see also Appendices 
O and El- It implies that the Noether energy-momentum tensor 

e'^^(X) = Tr (i)^^ J ^t;VF(X,p) + g^'^ (^-* - ^p-) + ^f/^^ ^ (23) 

is exactly conserved by the kinetic equation p2|) 

d^e^'iX) = 0. (24) 

Here potential energy density S^°^{X), which a probe particle with Wigner density F{X,p) 
would experience due to the interaction with all other particles in the system, is 

^pot(^) = Tr (i)^^ I ^ [ReE^F + ReG^E''^] . (25) 



Here in compliance with Eq. (O we define the factor (|)jj 



1/2 for neutral bosons (real fields) 
1 else 
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The interaction energy density £^'"*(X) specifies that part of the total energy density which 
is due to interactions. In simple cases it relates to S^°^ by a simple counting factor, namely, 
if all the interaction vertices of a theory have the same number rii of lines attached to them 

g'^^'^X) = -£P°*(X). (26) 

ni 

In particular, for two-body interactions one has = 4 and thus S™^ = while for the 

fermion-boson interaction rii = 3, which results in = |£^p°*. In Appendices O and iDl we 
discuss cases of this type. The additional term S^^^^^ ^ appears in Eq. (PHj) only in the case 
of derivative coupling, cf. Eq. ()B18|) . 

The considerations given above summarize the results of ref. Q| which are quite general. 
However they are restricted to local (point-like) interactions and void of derivative couplings. 
This excludes two important cases relevant to many areas in physics, nuclear physics in par- 
ticular. These are the cases of interaction mediated by finite range non-relativistic potentials 
and of derivative couplings like the P-wave pion-nucleon interaction. Since the considera- 
tions are rather technical they are exemplified in Appendices O and O There the results of 
[lfl| are generalized proving that also in these cases conserved currents and expressions for 
a conserved total energy and total momentum can be constructed. These two appendices, 
also provide further illustrations of the discussion given in the present section. 

The conserving feature are especially important for devising numerical simulation codes 
based on this kinetic equation. Indeed, if a test-particle method is used, one should be 
sure that the number of test particles is exactly conserved rather than approximately. For 
a direct application of this method, however, there is a particular problem with the KB 
kinetic equation. In the test-particle method the distribution functions are represented by 
an ensemble of test particles as follows 

F(X,p) ~ ^^'^ (X - ^^'^ - ' (27) 

i 

where the z-sum runs over test particles. Then the VF term in Eq. (fT^ just corresponds 
to the classical motion of these test particles subjected to forces inferred from ReS-^, while 
the collision term C gives stochastic change of test-particles' momenta, when their trajec- 
tories "cross". The additional term, i.e. the Poisson-bracket term |r™,ReG'^}, spoils this 
simplistic picture, since derivatives acting on the distribution function F appear here only 
indirectly. Namely they are encoded through derivatives of F™. This term is responsible 
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for backflow effects which restore the Noether current to be the conserved one. However, 
such backflow phenomena are difficuh to absorb into test particles, since they describe the 
response of the medium to the motion of the charges. In order to conserve the number of 
test particles between subsequent collisions, one would have to unite the additional term 
{r'^,ReG^} with the drift term VF even in the simplest case, when the collision term is 
charge conserving by itself (see Eq.()16|)) and derivative currents vanish, j^^gj,-) = 0. However, 
the interpretation of the additional term {F™, ReG^} causes problems within this picture, 
since it is not just proportional to the same ^-functions as in Eq. ()27p and thus can not be 
included in the coUisionless propagation of test particles. This problem, of course, does not 
prevent a direct solution of the kinetic equation. E.g., one can apply well developed lattice 
methods, which are, however, much more complicated and time-consuming as compared to 
the test-particles approach. 

Within the same approximation level the set of Dyson equations for Green functions 
G^^{X,p) provides us with an alternative equation for F 

MF - ReG^P'^ = ^ ({r,F} - {r°, A}) , (28) 

which is called the mass-shell equation, since in the quasiparticle limit it provides the mass- 
condition M = 0. This equation coincides with the kinetic one (1121) only within the first- 

nnnn 

order gradient approximation [Sjj lUl 112; 115|| ' while both equations are exactly identical before 
the gradient expansion. In view of this still remaining difference the practical recipe is to 
forget about the mass-shell equation since the retarded equation (fTUI) determines the 
spectral distribution, and to treat Eq. (fT^ as a proper quantum kinetic equation. Still 
this is an ambiguous recipe, which historically was one of the motivation to proceed to the 
Botermans-Malfiiet form of the quantum kinetic equation. 

D. Botermans Malfliet form 



As can be seen from the mass-shell equation (PHj) and Eq. ()12p 



BliilQ,li3 



the gain 



rate departs from FT /A only by corrections of first order in the gradients 

= rF/A + 0{dx) (29) 

(in equilibrium both equate each another). This fact permits to substitute the r.h.s. estimate 
for F™ in any of the gradient terms, i.e. in the Poisson-bracketed terms of Eqs. ()12|) and 
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and neglect the correction 0{dx) as it leads to terms of already second-order in the 
gradients. Upon this substitution, first proposed by Botermans and Malfliet Q|, one arrives 
at the following form of the kinetic and mass-shell equations 

F 



VF{X,p) 



r-,ReG^ 
A 



C{X,p), 



(30) 



(31) 



which are already exactly identical, as they were before the gradient expansion, and still 
equivalent to those in the KB form within the first-order gradient approximation. The so 
obtained equation (jHUj) will be called quantum kinetic equations in BM-choice. This equation 
exactly conserves the following effective current 



JBM-eff(^) = eTr J 



v''F{X,p) + ReJ: 



R 



OF 



ReG 



Rd{rF/A) 



(32) 



dp^ dp^ 

which differs from the Noether current in terms of the order of 0{dx), provided a $- 
derivable approximation is used for self-energies. All KB-choice properties of Eq. (j30|l 
within a ^-derivable approximation also transcribe to BM-choice through the substitution 
pin _ YF/A in Eqs. (PT|) . ()22|1 and (f^K|) . This substitution, however, touches the accuracy 
of those relations. For instance, the conservation laws of the Noether currents ()20|) and 
the energy-momentum tensor (PHj) are then only approximately fulfilled together with the 
corresponding consistency relations ()21|) and ()22|) which now look as 

d^p 



eTr 



/ WT' [^^'^'''^^ - {ReG^TF/A}+C] ^ d,j^,^^,, 



(33) 



p^d^p 



{Re E^,F} 



ReG^,^-^ \ +C 



TF 



n.b.. (2vr) 

respectively, and hold only up to first-order gradients. 
The effective BM-current (jH^ was used by S. Leupold 



A 



(34) 



121 as a basis for the construc- 



tion of a test-particle ansatz for the nonrelativistic case. In this case the additional term 
|rF/y4, Re G^} in the BM kinetic equation (j30|l is expressed in terms of the same distri- 
bution function as the drift term VF. Therefore, one can unify these terms to construct 
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equations of motions for test particles, which provide exact conservation of j^^-es- auto- 
matically fulfill this effective current conservation, the test-particle ansatz is introduced for 
the combination 

irA (l - ^ - F(X, p)~Y: (X - X.(T)) (p - p.(T)) , (36) 

rather than for the distribution function itself. Note that the energy p^{T) of the test particle 
is a free coordinate, not restricted my a mass-shell condition. W. Gassing and S. Juchem 

case. The equations of motion for 
the test particle, which follow from this ansatz, in particular give the time evolution of the 



off-shellness of a test particle 
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dT " r dT' ^"^^^ 

the origin of which can be traced back to the additional term jrF/A, ReG^} in the BM 
kinetic equation (jHUj) . Here M is the mass of the test particle relative to its on-shell value, 
see Eq. (|TT| . and this equation of motion implies that once the width drops in time the 
particles are driven towards the on-shell mass M = 0. This clarifies the meaning of the 
letter term for the off-shell BM transport: it provides the time evolution of the off-shellness. 

III. ENTROPY 

Another important feature of the kinetic description is the approach to thermal equilib- 
rium during evolution of a closed system. In terms of transport theory the sufficient (while 
not necessary!) condition of it is the existence of an H-theorem. Leaving aside all compli- 
cations associated with non-local effects in the collision term and possible lack of positive 
definiteness of the transition rates, discussed in ref. Isf, we confine our consideration to 
..pie app— io., cf. CD. A. de.on...ed . .f. 0. i„ the BM app.o..a«on to 
the quantum kinetic equations pUjl the H-theorem can indeed be formulated 

dAiX) = Tr 1 ^ In ^C"^''{X,p) > 0, (37) 

where the quantity 

Tr / — -Vt v^" n ■ 



'BM 



A 
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ReG^ I Tin 



F d 




(38) 



obtained from the l.h.s. of the BM kinetic equation ()30|) is interpreted as a entropy flow for 
the BM-choice. For the <l>-derivable approximation ()C1|) the r.h.s. of relation (j37|) takes the 
following form 

d'^Pi d'^p2 d'^Ps d^p4 



Tr 



d^n F 1 
-^ln-C^°^(X,p) =Tr- 



(27r)4 (27r)4 (27r)4 (27r) 



xi?i°^ (27r)^5^ (pi + P2 - P3 - P4) ( FiF2F3i^4 " F^F^F^F^ ) In 



F1F2F3F4 



(39) 



where is the transition rate determined by Eq. (jCMIl . This expression is indeed non- 
negative, since (x—y) ln{x/y) > for any positive x and y, and is of the second order in devia- 
tion from equilibrium {F-F^^), as both (F1F2F3F4-F1F2F3F4) and ln(FiF2F3F4/FiF2F3F4) 
approach zero at equilibrium. From the kinetic equation it follows that the deviation from 
equilibrium is of the first order in time gradients: {F — F^q) oc 0{dTF). This implies that 
the r.h.s. of relation (j37j) is of the second order in time gradients, which is, strictly speaking, 
beyond our first-order gradient approximation. However, from the point of view of practical 
use this feature is highly welcome as it guarantees equilibration. A further advantaged of 
the kinetic entropy flux (jnH|) is that in equilibrium its entropy merges that thermodynamic 
expression deduced from the thermodynamic potential in the $ derivable scheme jsl 3, Q| ■ 
In the case of the KB-choice fll2|) the situation is more controversial. Performing all the 
same manipulations with the KB kinetic equation |T2|l as those in ref. js], 
following relation 



we arrive at the 



d.sUX) = Tr 



where 



'KB 



Tr 



d^p 



(9ReS^ 



(40) 




- ReG^ Tin 



FdV' 
A dp. 



In 



Fdr' 
A dp^ 



(41) 



6H^ 



KB 



(42) 
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The KB entropy flow s^g is identical to the BM one Sgj^ up to zero-order gradients, while 
they differ in the flrst-order gradient corrections. One can easily obtain Sgp^ from the KB 
entropy flow by doing replacement in and a similar replacement in r°"*. 

The additional term 6H}<^b on the r.h.s. of relation is of the second order in gradients, 
due to the Poisson bracket and C'°'^ oc OldrF). Therefore, the r.h.s. of pn|) consists 
of two terms, which are of the same order of magnitude and one of them (5/7kb) is sign 
indeflnite. This prevents us from concluding the positive deflniteness of the r.h.s. of Eq. ()4Up . 
Alternatively we were not able to cast this term into a full divergence as to be included into 
the deflnition of the KB entropy flow. This fact by itself does not imply that the system 
does not approach equilibrium or even the absence of an H-theorem for the KB kinetic 
equation but suggests that equilibration should be tested in actual calculations. The local 
H-theorem we are looking for is a very stringent condition, providing monotonous approach 
to equilibrium. In fact, equilibration may well be nonmonotonous in time. 

Still, for the KB kinetic equation we are able to prove the H-theorem in a limiting case, 
i.e. close to local thermal equilibrium or for quasi-stationary state, which slowly evolve in 
space and time. To be deflnite, let us talk about the local thermal equilibrium. In terms of 
the distribution function 

F{X, p) = Fioceq. {X, p) + SF{X, p) (43) 

the above assumption implies that \5F\ <ti -Fioc.eq. and |c?x-^ioc.eq. I ~ Idx^F"]. Then we can 
write down 

^^jln^,ReG^ . (44) 



where 



Here, the remaining term 



j In -, ReG^ \ oc 0{5Fdx5F) + 0{5FdxF,,,,,^,), (46) 

can be neglected, as it has additional gradient smallness as compared to the flrst term 
on the r.h.s. of Eq. ()4()|1 . Here, we have taken into account that C'°'^ oc 6F and 
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|ln(F/F),ReG'^| oc (9x(Fioc.eq. + SF). Thus, from Eq. ^ we conclude that 

i^KB + "^-^kb) — ^^^^ local equilibrium, (47) 

which is the H-theorem for this case with the total entropy flow s^g + Ss'^^. Note that 5s^g 
is proportional to the collision term and hence equals zero in equilibrium. The applicability 
range of this result is the same as that for the memory entropy derived in for the BM- 
choice. 



IV. SUMMARY AND PERSPECTIVES 



In conclusion, we would like to summarize the present status of the two considered ap- 
proaches to off-shell transport. 

From a consistency point of view, the BM-choice looks more appealing, since it preserves 
the exact identity between the kinetic and mass-shell equation, a property inherent in the 
original KB equations Q. For the KB-choice this identity between the kinetic and mass-shell 
equation is only approximately preserved, namely within the validity range of the first-order 
gradient approximation. However, this disadvantage is not of great practical use, since, in 
any case, only one of these two equations, namely the kinetic one, should be used in actual 
calculations. 

For the construction of conservation laws related to global symmetries or energy and 
momentum the local collision term entirely drops out of the balance. Thus, the conservation 
laws solely depend on the properties of the first order gradient terms in the kinetic equation. 
In this respect the KB kinetic equation has a conceptual advantage as it leads to exact jl^ 
rather than approximate conservation laws, provided the scheme is based on ^-derivable 
approximations. Thereby the expectation values of the original operator expressions of con- 
served quantities (e.g., Noether currents) are exactly conserved. The reason is that the KB 
kinetic equation preserves certain contour symmetries among the various gradient terms, 
while they are violated for the BM-choice. Of course, within their range of applicability 
these two approaches are equivalent, because the BM kinetic equation conserves the charge 
and energy-momentum within the theoretical accuracy of the gradient approximation. Still, 
the fact that the KB-choice posses exact conservation laws put this version to the level of a 
generic equation, much like the Boltzmann or hydrodynamic equations, to be used as phe- 
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nomenological dynamical equations for practical applications. Such conserving dynamical 
schemes may be useful even though the applicability condition of the approximation might 
be violated at some stages of evolution. For instance, such a situation happens at the initial 
stage of heavy- ion collisions. As the conservations are exact, we can still use the gradient 
approximation, relying on a minor role of this rather short initial stage in the total evolution 
of a system. Moreover, exact conservation laws allow us to keep control of numerical codes. 

Although the KB kinetic equations posses exactly conserved Noether currents, a practical 
numerical approach (e.g., by a test-particle method) for its solution has not yet been estab- 
lished. The obstacle is the special Poisson-bracket term in the KB kinetic Eq. ()12|1 which 
lacks proper interpretation since the phase-space occupation function F{X,p) enters only 
indirectly through the gain-rate gradient terms. What is known is that this term encodes 
the backfiow component which ensures the Noether currents to be conserved. However such 
backflow features are difficult to be implemented into a test-particle scheme. This problem, 
of course, does not exclude solution of the KB kinetic equation, e.g within well adapted 
lattice methods, which are, however, much more complicated and time-consuming as com- 
pared to the test-particle approach. For the BM kinetic equation, on the other hand, an 
efficient test-particle method is already available 3, for the price that it deals with an 
alternative current rather than the Noether current. 

As a novel part we showed (cf. appendices) that the exact conservation laws in the KB 
kinetic equations, originally derived for local interaction terms which lead to a local energy 
momentum tensor, also do hold for derivative couplings and for interactions of finite range, 
like a non-relativistic potential. For the latter case, of course, only global conservation of 
energy and momentum can be achieved. In order to deal with the derivative coupling, we 
extended the <l>-derivable approach to this case. 

An important feature of kinetic descriptions is the approach to thermal equilibrium dur- 
ing evolution of a closed system. A sufficient (while not necessary!) condition is provided 
by an H-theorem. As was demonstrated in ref. at least, within simplest ^-derivable 
approximations for the kinetic equation in BM-choice an H-theorem indeed can be proven. 
The so derived kinetic entropy merged the equilibrium expression which in the context of 
derivable approximations results from the thermodynamic potential, cf. 0, 3, 3- -^^^ 
KB kinetic equation the result is by far weaker. Here we were able to prove the H-theorem 
only within simplest ^-derivable approximations and for a system very close to almost spa- 
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tially homogeneous thermal local equilibrium or stationary state. These results, in general, 
do not imply that the system does not approach equilibrium but suggests that equilibra- 
tion should be tested in actual calculations. Furthermore, the local H-theorem with a local 
entropy current which we considered for the BM-case may be by far too restrictive, provid- 
ing monotonous approaching to equilibrium. In fact, for kinetic equations with memory or 
non-local effects equilibration may well be nonmonotonous in time. 

Though the discussion in this paper is confined to problems of ^-derivable off-shell trans- 
port based on the first-order gradient expansion, essential progress has recently been achieved 
also in solving KB equations directly without any gradient expansion for selected examples. 
These concern non-equilibrium processes in scalar and spinor-scalar models on 1 and 3 space 



dimensions, see e.g. 
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2ll | that after a 



I and references therein. It was found 
comparably short but violent non-equilibrium evolution the time dependence of the Wigner 
transformed spectral function becomes rather weak even for moderate coupling constants. 
During this slow evolution the system is still far away from equilibrium. This fact provides 
a necessary condition for a successful gradient expansion and hence indicates a wide range 
of applicability of the approaches discussed in this paper. Even though the rapid far-from- 
equilibrium dynamics is formally beyond the scope of applicability of the gradient-expanded 
quantum kinetics, nevertheless, the KB choice includes all the ingredients required for such 
a treatment, i.e. the proper mean-field dynamics, together with the off-shell transport of 
particles, thereby satisfying exact rather than approximate conservation laws even far away 
from equilibrium. 

Further progress in understanding the properties of ^-derivable approximations to finite- 
temperature quantum field theory was reported concerning the question of renormalizabilty. 
The new results are equally applicable to quantum kinetic equations, both in KB or BM 
choices. In ref. it was shown that truncated non-perturbative self-consistent Dyson 
resummation schemes can be renormalized with local counter terms defined at the vacuum 
level. The requirements are that the underlying theory is renormalizable and that the self- 
consistent scheme follows Baym's ^-derivable concept. This result proves that there is no 
arbitrariness in studying the in-medium modifications of model parameters like the mass and 
the coupling constants within this class of approximation schemes. It is sufficient to adjust 
them in the vacuum, for instance, by fitting them to scattering data, in order to predict 
their changes in medium without ambiguity. This result also guarantees the standard 
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derivable properties like thermodynamic consistency and exact conservation laws also for 
the renormalized approximation schemes to hold. In ref. j^l the theoretical concepts 
for the renormalization, devised in [2^, were applied to the 0^-model, demonstrating the 
practicability of the method. 

In general, the symmetries of the classical action which lead to Ward-Takahashi identities 
for the proper vertex-functions are violated for ^-derivable approximations for functions 
beyond the one-point level, i.e., on the correlator level. This causes problems concerning the 
Nambu-Goldstone modes j^l in the broken symmetry case or concerning local symmetries 
(gauge symmetries) [2^1 on a level where the gauge fields are treated beyond the classical 
field approximation, i.e. on the propagator level. In ref. j3] it was shown that on top of any 
solution of a ^-derivable approximation which is constructed from a symmetric Lagrangian 
there exists a non-perturbative effective action which generates proper vertex functions in 
the same sense as the IPI effective action. These external vertex functions fulfill the Ward- 
Takahashi identities of the underlying symmetry. However, in general they coincide with the 
self-consistent ones only up to one-point order. Thus, usually the so generated external self- 
energy and higher vertex functions are different from the <l>-derivable expressions. Therefore, 
the pleasant property of the ^-derivable approximations, namely the conserving one, proves 
to be lost. The derivation of approximation schemes that fulfill all symmetry properties of 
the underlying classical action and at the same time are fully self-consistent and conserving 
still remains as an open task. 

As has been already mentioned, the gauge invariance may be lost in ^-derivable ap- 
proximations, too. In particular, this problem prevents applications of <l>-derivable approx- 
imations (including kinetic ones) to description quark-gluon plasma, based on QCD. This 
occurs because, in general, solution for dressed propagators and vertices do not satisfy Ward- 
Takahashi identities. This pathology shows up as an explicit dependence of results on the 
choice of the gauge condition. In ref. j^l it was demonstrated, in fact, that ^-derivable 
approximations have a controlled gauge dependence, i.e. the gauge dependent terms appear 
at orders higher than the truncation order. Furthermore, using the stationary point ob- 
tained for the approximation to evaluate the complete 2PI effective action boosts the order, 
at which the gauge dependent terms appear, to twice the order of truncation. This is still 
not a solution of the gauge problem in the rigorous sense but certain progress to its better 
control and understanding. 
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Appendix A: MATRIX NOTATION 

In calculations that apply the Wigner transformations, it is necessary to decompose 
the full contour into its two branches — the time- ordered and anti-time- ordered branches. 
One then has to distinguish between the physical space-time coordinates x, . . . and the 

corresponding contour coordinates x'^ which for a given x take two values x^ = {x^ ) and 
x~^ = (x+) {/J, e {0, 1, 2, 3}) on the two branches of the contour (see figure 1). 



t 



to t+ oo 

Figure 1: Closed real-time contour with two external points x,y on the contour. 

Closed real-time contour integrations can then be decomposed as 

f poo rto poo poo 

dx . . . — dx . . . -\- dx . . . — dx . . . — dx . . . , 

Jc Jto Joo Jto Jto 



(Al) 



where only the time limits are explicitly given. The extra minus sign of the anti-time-ordered 
branch can conveniently be formulated by a { — h} "metric" with the metric tensor in { — h} 
indices 

= K) = (J_°J (A2) 

which provides a proper matrix algebra for multi-point functions on the contour with "co"- 
and "contra" -contour values. Thus, for any two-point function F, the contour values are 

defined as 



F'^{x,y) := F{x\y^), i,je {-,+}, with 
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Fj{x,y) := F%{x,y) := F'\x,y)akj 



Fij{x,y) := a,k(rjiF^\x,y), = 5, 



ik 



(A3) 



on the different branches of the contour. Here summation over repeated indices is imphed. 
Then contour folding of contour two-point functions, e.g., in Dyson equations, simply be- 
comes 



H{x\y') = W'{x,y) = [F®G] 
= J dzF{x\z)G{z,y'' 



ik 



dzF\{x,z)G^\z,y) 



(A4) 



in the matrix notation. 

For any multi-point function the external point Xmax, which has the largest physical time, 
can be placed on either branch of the contour without changing the value, since the contour- 
time evolution from to x^^^ provides unity. Therefore, one-point functions have the 
same value on both sides on the contour. 

Due to the change of operator ordering, genuine multi-point functions are, in general, 
discontinuous, when ever two contour coordinates become identical. In particular, two-point 
functions like iF{x,y) = (TcA{x) B {y) \ become^ 



iF{x,y) 



^ iF~'{x,y) iF-+{x,y)^ 
yiF+'{x,y) iF++(x,y) j 



I (TA{x)B{y)) T{B{y)A{x)) ^ 



\ 



A{x)B{y) 



T-'A{x)B{y)) ) 



(A5) 



where T and are the usual time and anti-time ordering operators. Since there are 
altogether only two possible orderings of the two operators, in fact given by the Wightman 
functions F~^ and F^~ , which are both continuous, not all four components of F are 
independent. Eq. ()A5|) implies the following relations between nonequilibrium and usual 
retarded and advanced functions 



F\x,y) 



F\x,y) 



F-{x, y) - F-+(a;, y) = F+-{x, y) - y) 
e(xo-yo) {F+'{x,y)-F-+{x,y)), 
F-{x, y) - F^-ix, y) = F'+ix, y) - y) 
-Qiyo-xo) {F+-{x,y) - F-+{x,y)) , 



(A6) 



Frequently used alternative notation is F< = F + and — 
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where Q{xq — yo) is the step function of the time difference. The rules for the co-contour 

functions F etc. follow from Eq. ()A3|) . 

For such two point functions complex conjugation implies 

{iF-+{x,y)Y = iF-+{y,x) iF'+(X,p) = real, 

{iF+-{x,y)y = iF+~{y,x) iF+-{X,p) = real, 

{iF-{x,y)y = iF++(y,x) (iF-(X,p))* = iF++(X,p), 

{F^ix,y)y = F^iy,x) ^ (F^(X,p))* = F^(X,p), (A7) 

where the right parts specify the corresponding properties in the Wigner representation. 
Diagrammatically these rules imply the simultaneous swapping of all + vertices into — 
vertices and vice versa together with reversing the line arrow-sense of all propagator lines in 
the diagram. 

Contrary to the common case ()A5|) . the symmetrized contour convolution 

E{x)= j dz[F{x,z)G{z,x) +G{x,z)F{z,x)] (A8) 

is continuous, when two contour coordinates become identical. This can be easily checked, 
proceeding from relations ()A6|) for F and G functions. Moreover, for this symmetrized 
convolution with two coincident points we obtain a very simple expression in the Wigner 
representation, if all gradient corrections to the convolution are neglected (so-called local 
approximation), 

E'^'^X) = / [F-{X,p)G-{X,p) - F(X,p)++G++(X,p)] . (A9) 

In particular, this form is applicable to the potential (j25p and derivative ()D4|) energy den- 
sities. 

Appendix B: DERIVATIVE COUPLING 

To be specific, we consider a multicomponent system with different constituents a de- 
scribed by nonrelativistic fermionic and relativistic scalar bosonic field operators, summa- 
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rized as = {(f)a{x)}. The free Lagrangians of these fields are 

i (i4)ldt4>a - \dt4>l4>a ^V0lV0a ) nonrel. fermions 

\ (^9^(1) a ■ d^<Pa — ""^a^aj neutral rel. bosons (Bl) 



dfjL<p\d^ <pa — ^a4'i4'a charged rel. bosons 



We assume that these fields interact via linear derivative coupling, such that the in- 
teraction Lagrangian does not only depend on these fields but also on their derivatives 
£mt _ £int|0^^ (j)l^,d'^(f)a,d^(f)\}. The variational principle of stationary action determines 
Euler-Lagrange equations of motion for the field operators 

A* — 7 \ 7 \~ ~ — 7 \ — 7 \" ~' — ~i V^^) 



and the corresponding adjoint equation, where the "variational" 5-derivative is defined as 

Jf(^)---'-=^)---~^^{wjW)''')' ^^^^ 

This is the key definition, which allows to recast the local-coupling formulas to the derivative 
coupling case. In fact, the "variational" 5-derivative specifies the /it// derivative over /(x), 
implying that all derivatives acting on f{x) in the action should be redirected to other terms 
by means of partial integration before taking variational derivatives of f{x). 
The equations of motion can also be written in the differential form 

0(l)!i{x) 

and similarly for the corresponding adjoint equation. The J ai^) operator is the local source 
current of field a with mass m^, and Gq^{x) is the free evolution operator^ 

-dnd^ — ml for rel. bosons 
Go\x) ={ 1 (B5) 

idt H d for nonrel. particles 

2m„ 

with free propagator Go{y, resolvent, cf. Eq. 



^ Note that the first line in (|B5|) is not the nonrelativistic hmit of the second one. We have already divided 
the second line by 2ma, to take into account different normalizations of relativistic and nonrelativistic 
wave functions. 
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Invariances of the Lagrangian provide a set of conservation laws, the most prominent of 
which are those for the energy-momentum and certain currents. In addition to the standard 
canonical ener gy-m oment um tensor different representations of this tensor have been 
considered j29l . l3fll |. Using the Euler-Lagrange equations of motion and the definition of 
the source current ()B4|1 . one can show that the following form also defines a conserving 
energy-momentum tensor 



1 




x=y 



(B6) 



Here we have introduced the operators of the interaction-energy density £ and the 
potential-energy density 8^°^ 



(B7) 



^'^°'(^) = -\Y. (i)n.b. (^l(^)^a(x) + J, 



m 



Furthermore, we have singled out the contribution 



f) r'lnt f) r int 

0[da(t)a 




(B9) 



arising in the case of derivative coupling. Here and below the case of neutral bosons results 
from equating 0^ = in all the formulas. Proper counting is assured by the extra (D^^j^ 
factor which takes the value 1/2 for neutral boson (real fields) and 1 for complex fields. 

If the Lagrangian is invariant under some global transformation of charged fields (with 
the charges e^), e.g.. 



(BIO) 



there exists a Noether current defined as 



1 At 



dC 



d[d^(l)a 




J (conv.) 



(der.)- 



(Bll) 
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which is conserved, i.e. = 0. Here, we have decomposed it into two terms: the 

conventional one 

which is associated with the free Lagrangian, and the derivative term 

which is non-zero only for derivative coupling. 

To define the ^-functional for the case under consideration, all the steps described in 
ref. jr^ should be repeated. Then we arrive to the ^-functional that depends also on the 
gradients of mean fields {d^cpa and d^cpa) and Green functions {di^G{x,y) and d^G{x,y)) 
rather than on their values only. The variational rules of this functional formally look 
similarly to those in ref. jl^l 

d(pa [X) 

2 for real fields 

-i^a{x,y) = (T) T^TTTT-T X { (B15) 
oi'^ayyjX) I I for complex fields 

but should be understood in terms of the variational J-derivative of Eq. ()B3|) for one-point 
functions (like (f)a{x) and A(x)) and its generalization 

5iG{y,x)- 5,iG{y,x) ^ \5,{dtiiG{y,x)) ) ^ \5,{d^iG{y,x)) 

to two-point functions. Here 5q/ 5f)iG{y,x) means the conventional variation over G{y,x), 
which does not touch dl^iG, d^iG and dijidyiG terms in i$. Similarly to the variational 5- 
derivative of Eq. ()B3|1 . the 5- variation of Eq. ()B17|) means the full variation over G{y,x), 
implying that all derivatives acting on G{y, x) in the $ functional should be redirected to 
other terms by means of partial integration before taking variation over G{y, x). The factor 
A(x), appeared in Eq. ()B16|) . is an auxiliary scahng parameter of the coupling constant. 
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In terms of the $ functional the additional derivative contributions to mean values of the 
energy-momentum tensor ()B9|) and current ()B13|) take the following form 

5$ „ 5$ 



*^(der.) 



+ 



^ (der.) 



( (2)n.b. 



5$ 



5(9-iG,(z,x)) 



■ d^AGaiz, x) + d'liGaix, z 



5$ 



5(a-iG,(x,^)) 



(B18) 



'(der.) 



(dor.) 



-i^e. 

a 



5$ 



5$ 



+ 



d^ 



iG'a(a;, 



5$ 



(5$ 



5(6'-iG',(x,;2)) 5(6'-iG'„(^,x)) 



■ iGaf^jX) 



(B19) 



while the rest terms of 0^*^ and retain the same form as that for local coupling, cf. ref. 



17[. Here the variation is also understood in terms of Eq. ()B17|) to take account of the 
d^dyiGa{y,x) dependence of the $ functional. 

The next step to the kinetic description consists in gradient expansion of Kadanoff-Baym 
equations and all the related quantities. Expansion of the equations of motion up to the first 
order in gradients implies that the conserving quantities and self-energies, except for possible 
memory terms in the collision integral, are required only up to zero order in gradients. 
These zero-order quantities are determined by the local $ functional, where all gradient 
corrections are neglected. Since in the local approximation d^'iG, dy'iG and d^dy'iG transform 
into — ig'^iG(X, g), iq'^iG{X^q) and —iq^iq^iG{X,q), respectively, no partial integrations are 
needed for the variations of Eqs. ()B15|l - ()B16j) . This means that conventional variation 
rules of Eq. (jH} still hold in this case. At the same time, derivative contributions to the 
conserving quantities, Eqs. ()B18|) and ()B19|) . involve only variations over derivatives of the 
Green functions and, hence, should be carefully defined within the local approximation for 
the particular application considered. 



Appendix C: NONRELATIVISTIC NUCLEAR MATTER 

Currently, calculations of ground-state and low-temperature properties of nuclear matter 
are performed within the G- or T-matrix approximations to the self-energy 0, Q, Q- 
Based on realistic nonrelativistic nucleon-nucleon potentials, they quantitatively reproduce 
phenomenological properties of nuclear matter. However, already for the ground state the 
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resulting chemical potential, i.e. the single particle separation energy, deviates from the 
binding energy per particle violating the Hugenholtz-van-Hove theorem. This is a mani- 
festation of problems with the thermodynamic consistency in these approximations, which 
gets even worse at nonzero temperatures. This problem was discussed in refs. 0, Q, 0| . 
A consistent way to overcome this problem consists in using a self-consistent T-matrix ap- 
proximation [31] based on the ^-derivable approximation. 

Dynamic simulations of the nuclear matter are even more demanding to the choice of 
approximation to the self-energy, because the requirement of charge and energy-momentum 
conservations should be met except for that of the thermodynamic consistency. Again 
all these requirements are met provided the approximation is <l>-derivable. Since dynamic 
simulations are much more complicated as compared to static ones, up to now they were 
performed in a simpler approximation to the self-energy, i.e. the direct Born approximation 
P, 0|, which provide qualitative description of the dynamics. These simulations were based 
on the Kadanoff-Baym equations without any gradient expansion. Here we would like to 
call attention to the fact that the use of the gradient expansion in the KB form (see subsect. 
Ill would simplify these dynamic simulations and at the same time preserve the pleasant 
features of exact conservations and thermodynamic consistency. 

In view of reasonable level of complexity feasible for current computing, we confine our 
consideration to the full Born approximation to the $ functional 



i$ 



HFB 




s/ 

111 



+ 




Born 



(CI) 



which includes the Hartree-Fock contribution (the first two diagrams in Eq. ()C1|) ) and 
the true Born contribution (the last two diagrams). Here the wave line symbolizes a 

nonlocal nucleon-nucleon potential V{\ xi — X2 |), or y(| q |) in the momentum represen- 
tation. For simplicity, below we denote the latter as V{q), keeping in mind that in fact it 
does not depend on neither go nor direction of q. 
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Note that the (^^"'"^ part gives rise to the self-energy containing internal vertices. This 
implies that the corresponding collision term involves non-local effects (see discussion in ref. 

ri 

However, only "spatial non- locality" appears in the collision term, while the memory 
in time is absent since V{\ Xi — X2 \)S{ti — ^2) is time local. According to the general 
consideration of ref. exact conservations in the gradient approximation take place if all 
the non-local terms are consistently taken into account up to first-order gradients. Below 
we show that in the particular case of the <l>^^^ functional the exact conservations hold true 
even if we neglect the "spatial non-locality" generated by <|>^^^. These exact conservations 
imply the global rather than local conservation of the energy-momentum, which is in fact 
natural for the case of instant interaction of finite range considered here. 

Neglecting the gradient terms induced by the finite range of V, we consider the 
functional in the local approximation, where all Green functions in the Wigner representation 
are taken at the same centroid coordinate X. Alongside in some variational expressions we 
use an X-dependent local $ functional, $(X), where the last spatial integration is omitted, 
i.e. 

$ = y" dX$(X) (C2) 
The <l>^^^ of Eq. gives rise to the following local collision term 

^HFB-loc _ f dV2 d% d% HFB 

J (27r)4 (2vr)4 (2vr)4 
X [FiF^F^F^- FiF2hF^5\pi+p2-pz-pA), (C3) 

^HFB ^ (2^ ^y^^^ _ ^ y^p^ _ ^^)]2 ^ (^4^ 

where Fi = F(X,pi), etc., cf. Eq. (fTTI) . 



1. Charge Conservation 

In the G and G^^ Green functions are encountered only in different H — 

diagrams, and hence we can vary G and G'^^ independently. Therefore, is invariant 

under the following transformation 

G-(X,p) ^ G-(X,p + e(X)), G++(X,p) ^ G++(X,p-e(X)) (G5) 
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with F, F and V kept unchanged. Here E,{X) is an arbitrary function. If |^(X)| <^ 1, 
transformation ()C5|) reads 

5G-~=e,(X)^, 5G^^ = -UX)^—. (C6) 

Performing variation of c^^™'^ under the transformation ()C6|) within the canonical variation 
rules 0, we arrive at 

Here we have used that the integral 

dV f^dReG^ ^ 1 ^ /• d^ ^ndG^ 

11 I I 1 „ \ J t I — I -1 . , I -»> I \ *■ 



Tr / r ; + ReS^^ = — Im Tr / — ^-E^^^ = 0. (C8) 

J (27r)4 V dp, dp J 2 J (27r)4 ^ > 

equals zero due to analyticity of G^ and T^^. Thus, we obtain the relation 

J {2T,y \ dp„ dpj ' ' ' 

which guarantees the Noether current conservation, cf. Eq. ()2ip with j^^gj,-) = 0. 



2. Energy Momentum Conservation 

In order to construct the conservation laws related to space-time homogeneity we perform 
the following transformation 

G--iX,p) G-iX + aX),p), (CIO) 

with F, F and V kept unchanged. This transformation only acts on <|)^^^~ and <|)^^^+, i.e. 
those involving only " — " or " +" vertices, respectively, 

5$HFB _ J ^x^,{X)d'' - <I>H™(X+)) (Cll) 

where are understood in the sense of ()(^2j) . Note that 

i - = iTr / T^S^^^F 

J {2n)'^ 

+ ^iTr J ^(™eG^ + (ReS«-E"» 

- Tr y ^ (FRe + (Re - E^^) A) , (C12) 



29 



where the last integral is again zero due to analyticity, similarly to ()C8|) . Here the first 
term on the r.h.s. results from the first two (Hartree-Fock) diagrams in Eq. ()C1|) . while the 
last two integrals follow from the last two (Born) diagrams. Alternatively, we can perform 
variation of applying the canonical variation rules (jH)) 



m 



HFB 



dXe^(X)Tr 
2i / dXUX)TT 



d^p 
(2^ 



iS__ 



diG- 



-iS_ 



dF 



dX„ 



dX„ 



-i / dXUX)TT 



=0 



where the last integral is again zero due to analyticity. Therefore, we arrive at the important 
identity 



Tr 



d^p 



dX„ 



dX„ 



(2vr) 



(C14) 



Next we investigate the following transformation 



G^-{X,p) ^ G-(X,A^,(XK), G++(X,p) G++{X,A-J{X)p^) (C15) 



(C16) 



for the entire while F, F are kept unchanged. For the Hartree-Fock part <l>^^ one 

finds 



dV dV = (detA)-' d^pdV 



y dX [(detA)-'-l] <I>^^(X- 
j dX [(detA)2- 1] $«^(X+), 



(C17) 
(C18) 

(C19) 



where again <|>^^~ and are the diagrams involving only "— " or "+" vertices, 

respectively. In general, an arbitrary diagram <|)"g."a{-)^ consisting of no Green functions, 
nx of " — " interactions and no " +" interactions, transforms as 



§^nG,n^{-) = J dX [(det A)-("«-"^+^) - 1] $"g/"a(-)(x), 



(C20) 
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since the change of each momentum integration gives (det A)~^ (cf. Eq. ()C17|) ). while the 
transformation of the 5 (p) -function at each vertex (apart from one vertex due to global 
momentum conservation of the diagram) produces det A 

5\p + p' -q) =detK5\p + p' -q). (C21) 

Similarly, 

5$nG,n,(+) = j dX [(det A)"«-"^+^ - 1] <|."g.«a(+)(x). (C22) 

According to these rules the "— " and "+" diagrams of the second order in the interaction, 
^Born- g^j^j (|,Born+^ transformed as follows 

j dX [(det A)-3 - 1] (X), (C23) 

j dX [(det A)=^ - 1] <I>^°™++(X). (C24) 

If the A-transformation is infinitesimal, A^j,(X) = 1 + uj^y{X) with | uj^^, 1 and det A = 
1 + Trci;, det(A^"'^) = 1 — Trcj, we obtain 

8 ($HFB- ^ ^HFB+^) ^ _ j ^X2TlUJ (<I>^^-(X) - $«^+(X)) 

- j dXSTrcu (<I>^°™-(X) - <I>^°™+(X)) , (C25) 

cf. Eq. ()C12|) . The $Born-+ ^BoTn+- components, i.e. those containing both "-" and 
"+" vertices, are modified by only the V transformation. Moreover, this transformation 
leaves them invariant 

J J (27r)4^ V '^i^- dp^ S^V+ dp, J 

f ^ f dV , /5i$^°™"+(X)^, 5i$i^°™-+(X)^ A 1 dV ^ 

since ^^^'^ ^ is symmetric with respect to V~ and . Thus, 

5i$«™ = - j dX(Tra;)iTr j ^ [2^^^^ 

- ^ (r^'Re + (Re - . (C27) 

Alternatively, we can perform variation of applying the canonical variation rules 
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2i / dXu 



Tr / S^p^ ( r-^^ + Re E^|^ ) + Q-(X) 



dp. 



dp. 



(2vr)^ 



dp. 



dp^J' 



(C28) 



=0 



where the last integral again equals zero due to analyticity. Here we have introduced the 
quantity 



2ig^'^(X) = Tr / fP^y 5i$«™+(X) diV 



(27r) 



6iV- 



dpu 



6iV+ 



dp. 



(C29) 



arising from variation over V. All we have to know about this quantity is that it is Q'^'^ = 
when /i = and/or u = 0. This property results from po-independence of V{\ p |). In 
particular this property yields 



J d=^x9^g^"^(x) = 0, 



(C30) 



where the X-integration runs only over space. 

Hence, comparing Eq. ()C28|) to Eq. ()C27|) . we arrive at another important identity 



Tr 



( pn^^ ^ |F ■ J ^ ^^^^^^ 



(2vr) 



dp. 



d'^p 



4 



^P^^ 
d^p 



/ 7^ (r°ReG^ + (ReS« - E"'')F) . (C31) 
J (27r)* 



(27r)4 4 y (27r)4 

We turn now to the r.h.s. of the consistency relation for energy-momentum conservation 



(1221) 



(C32) 



integrated over space, which is suitable for the global conservation. In this expression the 
local collision term ()C3|) drops out according to Eq. (fT^ . It can be transformed by means 
of the identity 

d'^p 



-y {v. /} 



where ip and / are arbitrary functions, with the result 



(C33) 



Tr / d^X^'^ 



(27r 



-Tr / d=^X-^ (Re E^d'^F + T'^'d^Re G^) . 
J (27r)4 



(C34) 



32 



Now, applying identities (jnil) . and KJM^ to the r.h.s. of Eq. ([HHill . we obtain 

= J d'^Xd^g^^'Ti J (TOeG^ + ReS^F) (C35) 

which is precisely needed for the global conservation of the Noether energy-momentum 

d^Xe°^(X) = 0, (C36) 
since for the case under consideration S^^^ = cf. Eq. (j26p . 



d_ 



Appendix D: NUCLEON-PION SYSTEM 

For the discussion of the physical aspects of the nucleon-pion problem we refer to ref. 



j^. Here we would like to clari 
of pion-nucleon interaction 



y some technical details. We choose the nonrelativistic form 



2'-' = g^^ [(o-- V)(t- 0)] ^, (Dl) 

where 'ip and are non-relativistic nucleon and Klein-Gordon pion field operators, respec- 
tively. We accept a simple approximation defined by the following $ functional 

N 



2 T T 2 J (27r)4 (27r)4 (27r)4 ^' 
N 

X iG^^(X,pi) {-ig){-iq-a)TiA,,{X,q) {-ig){iq ■ <t)t iG^\X,p,), (D2) 

where G and A are the nucleon and pion Green functions, respectively, and Tr runs over 
spin and isospin indices. Here we have assumed an isotopically symmetric system, where 
the pion Green functions of all isotopic charges coincide. Though this approximation is 
evidently oversimplified to produce quantitative results, cf. ref. jsj, it is able to reproduce 
qualitative features of the dynamics. Moreover, this approximation is at the edge of the 
present computing abilities. The formal basics of the ^-functional formalism are given in 
Append. IHI 

The charge current, defined by Eqs. (j2U|) and (jBllI) . relates to the baryon number 
conservation and hence is trivial from the point of view of the pion-nucleon interaction. 
Indeed, to prove the baryon number conservation we should perform transformation ()G5|) 
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with F]\f, Fjy and A*-' kept unchanged. The pion Green functions Aij are not subjected 
to this transformation, since pions are neutral from the point of view of baryonic charge. 
All the subsequent considerations are completely identical to that of the Fock diagram (the 
second term in Eq. (jC^lj) ) and lead to the same final result ()(]9j) . i.e. to the exact Noether 
current conservation. 

The energy-momentum conservation is more instructive in this respect. Before proceed- 
ing to the conservation laws themselves we should define the derivative contribution to the 
energy-momentum tensor (|B18|) . In our case of vanishing mean fields the pion Green func- 
tion enters the ^-functional only doubly differentiated. Therefore the expression ()B18jl takes 
the form 



where 5q is already the conventional variation. In the Wigner representation with due regard 
for Eq. (jA9|l it transforms into 



Contrary to usual convention, here we use Greek indices /i and u for the components of 
3-vectors in order to distinguish them from the "H — " summation indices. The potential 
energy density is still determined by the standard expression but in terms of the sum 
over nucleons and pions. Within the approximation of Eq. ()D2|) . S^°^ can be alternatively 
expressed as 



where ^'^^ (<|)'^^+) refers to the ^'^^ functional with removed integration over dX and all 
vertices being of "— " ("+") type. In view of relation ^Hjl . 




(D3) 




X [iG-{X,p,) {-q''a^)T iA-(X,g) {iq ■ (t)t iG'-(X,p2) 

- iG'++(X,pi) (-gV^)TiA++(X,g) (ig ■ o-)t iG++(X,p2)] • (D4) 



(D5) 



(D6) 
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1. Energy Momentum Conservation 



We briefly repeat the steps proving the exact energy-momentum conservation for the 
nonrelativistic nuclear matter (Append. 2|1 with the emphasis on the specifics of the 
derivative coupling. 

First the transformation of Eq. ()C10|) for the nucleon Green functions together with the 
corresponding transformation of the pion Green functions 



A-{X,p) A-{X + aX),p), A++(X,p) ^ A++(X-e(X),p) 



(D7) 



with and F^^ being kept unchanged has to be performed. This transformation is unaffected 
by the derivative coupling, and in the similar way as before we arrive at the identity 



Tr 



(27r)4 



(D8) 



Here S and 11 refer to nucleon and pion self-energies, respectively, and subscripts N or n 
correspondingly attribute a quantity to either nucleon or pion subsystems. The r.h.s. of this 
identity is written with due regards to Eq. ()D6|) . 

Let us now perform the transformation ()G15|) for both nucleon and pion Green functions, 
as well as q factors encountered in vertices of Then the variation of ^'^^ gives 



5$" 



N 



dXuj 



flU 



_2gf^u ($-A^-(x) - $'^^+(X)) + 2£, 



fll/ 

(der.) 



(D9) 



where the 8^^^^-^ results from the variation of q factors in vertices of <l>'^^. Alternatively, 
performing variation of accordingly to the canonical variation rules and equating 
the result to expression ()D9j) . we arrive at another identity 

d^p 



Tr 



(2vr) 



^ dpn dpn ) 2 V dpa dp ^ 



(der.)' 



M / J 



(DIO) 



The r.h.s. of the consistency relation for energy-momentum conservation (j22j) now reads 

p'^d'^p 



K" = Tr 



{2nY 



[{{ReJ:^,F^}-{ReG^,r%}) 
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+ - ({Ren^F4-{ReA^,C 



})]> 



(Dll) 



By means of identity ()C33|) it is transformed to the form 




(D12) 



Now, applying identities ()D8|1 and ()D10|1 to the r.h.s. of Eq. ()D12j) . we obtain 



(der.) 



(D13) 



which is precisely needed for the local conservation of the Noether energy-momentum. 
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